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*5 ' Abstract 

In this letter, we determine the K-distribution function for a gas in the presence of an external 

field of force described by a potential U(r). In the case of a dilute gas, we show that the K-power law 

distribution including the potential energy factor term can rigorously be deduced in the framework 

^>. ■ of kinetic theory with basis on the Vlasov equation. Such a result is significant as a preliminary 

to the discussion on the role of long range interactions in the Kaniadakis thermostatistics and the 

m 

underlying kinetic theory. 

o. 

OO PACS numbers: 51.10.+y; 05.20.-y; 05.90.+m 

o" 
> 

•i-H 

X 
S3 



'Electronic address: jmsilva@astro.iag.usp.br 



^Electronic address: limajas@astro.iag.usp.br 
-t-Electronic address: raimundosilva@uern.br 



I. INTRODUCTION 

Over the last two decades, a great deal of attention has been paid to the so-called power- 
law statistics, both from theoretical and observational viewpoints. The main motivation is 
the lack of a comprehensive and rigorous treatment including gravitational and Coulombian 
fields, or more generally, any long range interaction for which the assumed additivity of the 
entropy present in the standard approach is not valid (see, e.g., [l|, ]2j, y|, U, ]5|, [g, ]7[ and 
references therein). 

In this context, the Tsallis |l|, Id, E, \i, U U IJ and Kaniadakis [si, E, U, U, QJ, U, U, UM, UM 
power-law statistics are the most investigated frameworks. In the former case, the standard 
Boltzmann-Gibbs formalism is extended through a new analytic form for the entropy S q = 
/cb(1 — YliiPi) I '(Q ~ 1)' where k-Q is the standard Boltzmann constant, p, t is the probability 
of the ith microstate, and q is a parameter quantifying the degree of nonextensivity. This 
expression has been introduced in order to extend the applicability of statistical mechanics 
to system with long range interactions and has the standard Gibbs-Jaynes-Shannon entropy 
as a particular limiting case (q = 1). Even the so-called g-nonextensive kinetic theory 
has also been developed and applied for many different contexts ranging from plasmas 



,« 
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5] 



to gravitational systems p]. 

On the other hand, recent studies on the kinetic foundations of the K-statistics also leads 
to a power-law distribution function and a K-entropy which emerges naturally in the frame- 
work of the kinetic interaction principle (see, for instance, Ref. [8|, |3, LlOJ] ) . Several physical 
features of the so-called /^-distribution have also been theoretically investigated, among them: 
i) the self-consistent relativistic statistical theory [l0(] , ii) the framework of nonlinear kinetics 
[111 ] and iii) the H-theorem and generalization of the chaos molecular hypothesis 121 ] . Actu- 



ally, this K-framework leads to a class of one parameter deformed structures with interesting 
mathematical properties [13( . In particular, the canonical quantization of a classical system 
[14| , and the so-called Lesche stability have also been discussed in the K-framework [15| . Still 
more important, a consistent form for the entropy (linked with a two-parameter deforma- 
tions of logarithm function), which generalizes the Tsallis, Abe and Kaniadakis logarithm 
behaviours [la] have also been found. 

In the experimental front, there also exist some evidence closely related to the behavior 
predicted by the ^-statistics, namely, cosmic rays flux, rain events in meteorology |l3|], quark- 



gluon plasma 17|, kinetic models describing a gas of interacting atoms and photons 18], 
fracture propagation phenomena 19( , and income distribution 20] , and even the possibility 
of improved financial models has also been investigated 21]. 

Mathematically, the K-framework is based on ^-exponential and ^-logarithm functions 
which are defined by 

exp K (/) = (v / l + ^ 2 / 2 + ^/) 1/K , (1) 

M/) = (/" - r K )/^ (2) 

ln K (exp K (/)) = «p„(ln„(/)) = /. (3) 



The K-entropy associated with this K-framework is given by 

S K (f) = - fd 3 pf[- 



fK f— K 



2k 



(4) 



which recovers standard Boltzmann-Gibbs entropy, S K= o(f) = — J fin fd 3 p, in the limit 



K — ■> (see Ref. 



9j for details). 
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The so-called Kaniadakis entropy reads 

S K = -Jd 3 pf\n K f = -(\n K (f)). 
while the equilibrium velocity distribution can be written as [8j, [9|, llO, 



(5) 
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fo(v) = A K 



1 + K 2 



rav^ 



2k B T 



mv 



2k B T 



(6) 



In this expression k B is the Boltzmann constant while the k index denotes a continuous 
parameter associated to the gas entropy, and whose main effect at the level of the distribution 
function is to replace the standard Gaussian form by a power law. The quantity A K is a 
normalization constant fixed by the total number of particles in a given volume. As it should 
be expected, the above expression reduces to the standard Maxwellian distribution in the 
limit k = for which A = (m/27i-k B T) 3/2 . 

In this letter we explore how the potential energy term can rigorously be introduced in 
the K-distribution 22J. More precisely, we deduce an analytical expression for the equilib- 
rium distribution of a dilute gas under the action of a conservative force field with basis 
on the stationary solution of the collisionless Boltzmann equation. As we shall see, this 
result is significant as a preliminary to the discussion of long range interactions according 
to Kaniadakis thermostatistics and the underlying kinetic theory. 



II. VLASOV EQUATION AND THE BOLTZMANN FACTOR 



In this section we discuss briefly the standard case, i.e., the kinetic description of a classi- 
cal gas under stationary conditions and immersed in a conservative force field, F = —VU(r). 
Typical examples are a gas in the earth's gravitational field and ions in an external magnetic 
field 
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241 ]. This kind of problem is important on their own because it permits to under- 



stand how the molecular motion is modified by force-fields different from those exerted by 
the containing vessel or even by the other particles of the gas. As widely known, its distri- 
bution function differs from the Maxwellian velocity statistics trough an extra exponential 
factor involving the potential energy whose general form reads 



f(r,v) 



n 



m 



3/2 



exp 



\mv 2 + U(r) 



(7) 



27tA;bT/ \ ksT 

where m is the mass of the particles, T is the temperature and no is the particle number 
density in the absence of the external force field. In addition, since this distribution function 
is normalized, it is easy to see that the number density is given by 

U(r 



n(r) = n exp 



hnT 



(8) 



where the factor, exp[— L^(r) /&bT], which is responsible for the inhomogeneity of f(r,v), is 
usually called the Boltzmann factor. Expression ([7j) follows naturally from an integration of 
the collisionless Boltzmann's equation 

F df 



<V +V .^ 

dt dr m dv 



0. 



(9) 



when stationary conditions are adopted along with the assumption that the total distribution 
can be factored 

f(v,v) = f (v)x(r), (10) 

where fo(v) represents the Maxwell equilibrium distribution function, and x( r ) is a scalar 
function of r. As one may show, after a simple normalization, the resulting expression for 
x(r) is exactly the Boltzmann factor for the potential energy of the external field, namely: 



x( r ) 



cxp 



k B T 



(11) 



and combining this result with equation (TTU1) the Boltzmann stationary distribution (EJ) is 
readily obtained. 



III. VLASOV EQUATION AND KANIADAKIS KINETIC THEORY 

Let us now consider a spatially inhomogeneous dilute gas supposed in nonequilibrium 
stationary state at temperature T. The gas is immersed in a conservative external field in 
such a way that f(r, v)d 3 vd 3 r is the number of particles with velocity lying within a volume 
element d 3 v about v and positions lying within a volume element d 3 r around r. In principle, 
this distribution function must be determined from the K-type Boltzmann equation: 

v.VJ - -V r UV v f = C k (f) (12) 

m 

where C K denotes the k-collisional term. The left-hand-side of the above equation is just the 
total time derivative of the distribution function. Hence, the effects appearing from the k- 
approach can be explicitly incorporated only through the collisional term. In particular, this 
means that the Vlasov equation, or the stationary Boltzmann equation takes the standard 
form (see Eq. (7)) 

v.V r / - -V r U.V v f = 0. (13) 

m 

In order to introduce the K-statistics effects we first recall that the factorizability condition 
as given by ffTUl) is modified in this extended framework. This means that the standard 
starting assumption must be extended. In the spirit of the K-formalism, a consistent k- 
generalization of ( TTU1) is 



f(r,v) =A K exp K 



\n K [j-) ln K x(r) 



(14) 



where Jq denotes the K-velocity distribution and the constant normalization has been in- 
troduced for mathematical convenience, and the functions exp K (/), ln K (/), were previously 
defined by Eqs. (JTJ and ©. 

Note that in the limit k — > the above identity reproduces the usual properties of the 
exponential and logarithm functions. In addition, since exp K (ln K /) = ln re (exp K (/)) = /, 
the standard factored decomposition ([TO]) is readily recovered in the limit k — 0. In what 
follows, the properties of ^-exponential and K-log differentiation 

d\n K f ( r- 1 + /-("+ 1 ) ^ df 

~lbT ~ V 2 ) Tx" (i5j 

dexp K (f) = exp K (/) df_ 
dx " ^1 + K 2 f 2 dx' 

5 



will also be extensively used. In particular, for the total ^-distribution (j!4p . we obtain 



V7 -/ x exp K ln K / (w)H-ln K x(r) 

V r f(r,v) = ^ pr -V r ln K x(r)x 

exp«[ln K / (w) + ln K x(^)J 



1 + 



2k B T 



[l + ^a^/o^ + ln^xW) 2 ] 1 / 2 



(17) 



exp£[ln K /o(u) + ln„ x(r)] \ k B T 
k ( ln K x(r 



v ,, x = exp K [ln K / (t;) + ln K x(^)] f mv ^ 



1 + 



2/c s T 



(18) 



[l + K 2 (ln K /oM+ln K x(r)) 2 ]V2 

Now, substituting V r / and V„/ given above into the stationary Boltzmann equation 
I2|) . and performing the elementary calculations one obtains 



V r In x ■ dr 



k B T 



W(r) ■ dr 



the solution of which is 



x(r)=exp K (-^ + C 



(19) 



(20) 



where C is an arbitrary constant. 

Here, inserting (120]) into (TH1) . and integrating the result in the velocity space it follows 

that 



K exp K 



hi K (^/o) 






+ C 



d v = n(r). 



(21) 



Finally, by substituting the expression of /o and considering a region where U(r) = 0, one 
finds 



M / exp. 



771V 



2A; B T 



+ C d 3 u = n , 



(22) 



and from normalization condition, uq = J fo(v)d 3 v, it follows that the unique allowed value 
for the integration constant is C = 0. Consequently, f[2"Uj) becomes 

U(r) 



X{r) = exp K 
which is the ^-generalized Boltzmann factor. 



k B T 



(23) 




FIG. 1: The K-velocity distribution function f K (r, v) for typical values of the k parameter on the 
interval [-1,1] (2 = A K ~ ). In contrast to Tsallis statistics [l|, |2(, this power law distributions does 
not exhibit a thermal cutoff, that is, a maximum value allowed for the velocity of the particles. 
The K-distribution also presents a striking mathematical property, namely, f K {j]) = f- K {rf). 



Finally, by inserting this result into ffMj) . we obtain the complete ^-distribution function 
in the presence of an external field 

-1 1/k 



f K {r,v) = A K 




mv z 



2k B T k B T 



+ K 



mv 



m_ 

2k B T k B T 



A K exp K (-E/k B T), 



(24) 



where E is the total energy of the particles. It thus follows that a generalized «;-exp factor 
for Kaniadakis' thermostatistics can exactly be deduced if the standard approach is slightly 
modified. In Fig. 1 we show zf(E) for some selected values of the ^-parameter (where 
z = A K ~ l ). Different from Tsallis power-law functions which can become finite for some 



values of the q-parameter, the ^-distributions are not finite. In other words, the thermal 
cutoff on the velocity space is not present in such distributions regardless of the values of k. It 
is worth emphasizing, however, that this ^-distribution presents the following mathematical 
behavior, viz., f K (rj) = f- K {rj). 

IV. CONCLUDING REMARKS 

In the last few years, several applications of the K-power law velocity distribution have 
been done in many disparate branches of physics [8-21]. However, such investigations are 
usually related with the k- velocity distribution function as given by equation (6). On the 
other hand, many physical systems involve naturally the presence of a conservative force 
field as happens for example with ions in a magnetic field. Probably, the most popular 
problem of a gas in a force-field is the planetary atmosphere. In the standard simplified 
treatment, the temperature is uniform and the tree-dimensional motion occurs under the 
action of a constant gravitational field along the z-direction. To all this sort of problems, 
the extended ^-distribution deduced here with basis on the Vlasov equation, namely 



f K {r,v) = A K 



-| 1/k 

mv 2 U(r)\ ( mv 2 U(r) 

' '" ' ~2k B T ~~ k B T ) +K \2k B T ~ ' TjT 



can be applied, and, might prove to be of extreme wide usefulness. Note also that a giroscopic 
term may also be added to the above power law distribution. In a rotating frame with 
constant angular velocity, the whole effect is just to add a Coriolis term —l/2muj 2 R 2 to the 
potential energy, where R is the perpendicular distance from the axis of rotation. In the 
Newtonian framework, such a term simulates a change in the potential energy due to gravity. 
Finally, it is worth mentioning that the present consistency between Vlasov equation and 
Xaniadakis thermostatistics also is valid in the context of Tsallis nonextensive framework 



22]. 
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